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Abstract 

We derive an exact expression for the two-point correlation func- 
tion for quantum star graphs in the limit as the number of bonds tends 
to infinity. This turns out to be identical to the corresponding result 
for certain Seba billiards in the semiclassical limit. Reasons for this 
are discussed. The formula we derive is also shown to be equivalent to 
a series expansion for the form factor — the Fourier transform of the 
two-point correlation function — previously calculated using periodic 
orbit theory. 



1 Introduction 



The statistical distribution of quantum energy levels is a much studied topic. 
It has been conjectured that generic, classically integrable systems give rise to 
uncorrelated quantum spectra while the energy levels of generic classically 
chaotic systems have the same statistical properties as the eigenvalues of 
random matrices || . This has been confirmed by semiclassical theory |3]. [| . 
and in a large number of numerical studies, but classes of systems have also 
been found for which it is not true; these include geodesic motion on surfaces 
of constant negative curvature ||, and the cat maps ||. 

Quantum graphs [7|, Q are mathematical models introduced in order to 
explore the connection between the periodic orbits of a system and the sta- 
tistical properties of its energy levels. The trace formula, in which the level 
density is connected to a sum over periodic orbits, is exact for graphs, rather 
than a semiclassical approximation, and the orbits can be classified straight- 
forwardly. However, despite the fact that numerical computations have re- 
vealed good conformance of the spectral statistics of many quantum graphs 
to the predictions of Random Matrix Theory (RMT), few conclusive analyt- 
ical results have been obtained so far. This is due to the fact that although 
some individual finite graphs can be shown to reproduce certain features of 
RMT behaviour || |K], [1T[] , the full RMT results can only be recovered in 
a limit in which one is forced to consider larger and larger graphs, and this 
necessitates finding general, combinatorial asymptotic techniques for dealing 
with the (non-trivial) length degeneracies of the periodic orbits. 

One family of graphs in which this goal has been achieved are the star 
graphs |l2j (defined below and shown in Fig. [I]), but in this case the resulting 
spectral statistics are neither RMT nor Poissonian (i.e. those of random num- 
bers). It turns out, however, that it is not the first time that such statistics 
have arisen in the connection with the study of quantum chaos. Our purpose 
here is to demonstrate that the star graphs have exactly the same two point 
spectral correlations as a large class of quantum systems, which we will refer 
to as Seba billiards. 

The original Seba billiard, a rectangular quantum billiard perturbed by 



as an 



a point singularity (also illustrated in Fig. |T]), was introduced in [[U| 
example of a system whose classical counterpart is integrable (the singularity 
affects only a set of measure zero of the orbits) but which nonetheless exhibits 
properties of quantum chaos. This construction was later generalized to all 
integrable systems jnj perturbed in the same way. We will refer to any 



2 



system in this class Seba billiard. 

The energy levels of a Seba billiard can be found by solving an explicit 
equation which depends on the levels of the original unperturbed system and 
on the boundary conditions imposed at the singularity. This equation takes 
the general form 

= 1, (1) 
where £(z) is the meromorphic function 

« 2) = E M^ (2) 

the sum being suitably regularized to ensure convergence. Here {Ei} are the 
eigenvalues of the unperturbed system, T/> n (x ) is the value of the nth unper- 
turbed eigenfunction at the position xo of the singularity, and the coupling 
constant A parametrizes the boundary conditions [13|, |14|j . Assuming that 
{Ei} are given by a Poisson process, one can then calculate the associated 
spectral statistics, such as the joint level distribution, asymptotics of the level 



spacing distribution |L4| , and the two-point spectral correlation function [15 
The results show the presence of spectral correlations but are substantially 
different from the RMT forms. 

Here we apply the methods developed for Seba billiards in ([TJJ to calculate 
the two-point spectral correlation function for star graphs, starting from an 
expression which is analogous to (0). The formula obtained will be shown to 
be a resummation of the expansion computed from the periodic orbit sum 
in [T2\. Our main result will be that this correlation function is the same as 
that already found for Seba billiards in the case when |?/; n (xo)| 2 = constant 
(e.g. when the billiard is rectangular with periodic boundary conditions) and 
A — > oo. We finish with a discussion of reasons for this coincidence. 



2 Quantum star graphs 

Star graphs are metric graphs of the type shown on Fig. [I] with a Schrodinger 
equation 

-■^*i = * 2 * i) XjGiO,^}, (3) 



dx 3 
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(a) (b) 

Figure 1: A star graph with v edges (a) and a Seba billiard (b). 

defined on the bonds and boundary conditions, for example 

*,-(<>) = * t (0), (4) 

E^OHo, (5) 

5 J 

^i(Li) = 0, (6) 

specified on the vertices. Here Lj is the length of the j-th bond, j — 1 . . . v, 
and the real variable Xj varies from to Lj, with corresponding to the 
central vertex and Lj to the outer vertex. The lengths Lj are assumed to be 
incommensurate; see [12| for further details. We refer to positive values of 



the parameter k for which the system is solvable as eigenvalues of the 

quantum star graph. 

Denoting the ordered sequence of eigenvalues by {h}^, we define the 
spectral density by 

oo 

d(k)=Y,8{k-ki). (7) 



The statistic we shall mainly be concerned with is the two-point correlation 



function 



R 2 (x) = ^(d{k)d(k+^J-5(x), (8) 

where d = (d(k)) is the mean density, 5(x) is the Dirac ^-function, and the 
average ( • ) is either over k, or over the bond lengths Lj (we shall specify 
which in each particular context). R2{x) is an even function and hence so is 
its Fourier transform, 

POO 

K(t) = 1 + 23? / (R 2 (x) - l)e 2niXT dT, (9) 
Jo 

which is called the form factor. 

A complete series expansion of the v — *■ oo limit of K(t) in powers of r 
around r = was derived for the star graphs in using the trace formula 
and a classification of the periodic orbits: 



OO OO i A 

4 J 



K{t) = exp(-4r) + £ £ -C jM r M ^\ (10) 

i=2 M=0 3' 



where 

> (A + J - 

(M + j-1)! 



Cj.M = (-2)" £ ^^iV'^V - " - K), (11) 



iC=0 



with 



F '(^' = iITWW' (12) 



and 



AT N 



Fj(K,N) = '^2^2F 1 (k,n)F j -i(K-k,N-n). (13) 



fc=0 n=0 



Explicitly, 

4. o s 32 4 16 = 128 6 16 7 64 s , „, 
K{r) = e" 4r + 8r 3 - yr 4 + -r 5 - — r 6 + -r 7 + -r 8 + o(r 8 ). (14) 
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K(X) 0.6 




Figure 2: The sum of the first 30 terms in the expansion (|H]) (dashed line), 
which converges in the range r < r* ?a 0.63, compared to the results of 
a numerical computation || of K(r) (circles). Also shown are the Pade 
approximations to the series of order [21/20] (thin solid line) and [23/23] 
(thick solid line). 



In this calculation, the average in (§) was over k. The result is in excellent 
agreement with the numerical data (see Fig. ^) but is limited by the fact 
that the radius of convergence of the series is finite, being approximately 0.63 
(found by applying Cauchy's test to the coefficients in the series, but see also 
Fig. HD. The range of convergence can be extended using Pade approximation 
(again, see Fig. |2|), which suggests that the singularity causing the divergence 
is not on the positive real line [16| . 

Here we approach the problem from a different direction: it is possible to 
solve equations (|3])-(||) to derive an explicit condition on k to be an eigenvalue. 
Indeed, the general solution of (|3|) on a star graph can be written in the form 
\F,(x) = Aj cos(k(x + 4>j)), j — 1, . . . , v. Applying condition ([]), we obtain 
<pj = —Lj while condition (^) on the central vertex implies Aj cos(Ljk) = 
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const. Finally, applying condition (^j) and dividing by Aj cos(Ljk) we obtain 

V 

^2tanLjk = 0. (15) 
i=i 

Similar expressions can easily be found when different boundary conditions 
are applied at the central vertex. The general equation reads 

- 1 
^^tanLj/c = — , (16) 

i=i 

where A is arbitrary parameter. However, in the limit as v — > oo, A fixed, 
the two-point correlation function turns out to be independent of A (see 
the comment following equation ((49|)). Our calculations will therefore be 
performed for A" 1 = 0. 

Note the similarity between ([16]) and the quantization condition (p]) for 
Seba billiards when |^ n (x )| 2 = constant. 

Condition ([T^) means that k is an eigenvalue if and only if it is a zero of 
the function F{k) = ^J =1 tanLj/c, and so we can express the density d{k) 
as 

d(k) = — f \F\k)\e izF{k) dz = — [ V L 2 S - e iz ^U^ k dz. (17) 

27T / 27T / COS L ck 

Our analysis of the spectral correlations will be based on this representation. 



3 Mean density. 

As an example of the techniques to be employed later, we begin by calculating 
the mean density d defined as 



d= lim (d(k)) (18) 

AL->0,fc->oo \ / {Lj} 

where now the average is with respect to the individual lengths of the bonds, 
rather than over k: 

r L ^ L r L ^ L du dLv 
(• W = / ••• / • al" aT' (19) 



i) 
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That is, we assume that the lengths are independent random variables dis- 
tributed uniformly on the interval [L ,L + AL]. We also assume that AL 
and k tend to their respective limits in such a way that ALk — » oo. 
Applying this averaging to flT7|) we obtain 



1 P , f f L °+ AL e^EU^^i dLi dL v 



d(k) hL j} - ^y^^Ey -J L ° COS 2 kLs AL 



v r dz ( [ L » + ^ nkL dL\ f * dL 



2W-oc \Jl \J Lo cos 2 kLAL 

v 

2^ 



oo 

v-1 



r\z)~g{z) dz. (20) 



L +AL AztaxikL jf j /-tan k(L +AL) 

r e aL . . M) / iztwkL 



Here 

where we were able to approximate L by Lo because it is slowly varying 
(compared with tan kL) and ultimately we will take the limit AL — > 0. Now, 
since tan kL is a periodic function with the period of ir/k, and the integration 
is performed over the interval containing approximately ALkjix periods, we 
can further approximate 

g( z ) = ^(— T e itM dtfmkL + 0{l)) 
ALk V. 7t J 

« 2L 5(z), (22) 
where 0(1) is a quantity which is bounded as kAL — > oo. Similarly, 

L ° +AL iztankL^L Lq f^^^^jztankL d tan kL 



/(*) = e — — = XTT / 



AL ALA; 7 tanfcLo 1 + tan 2 fcZ, 



1 



oo 



——da = e-\ z \ (23) 
7T J_oo 1 + a 2 

where the last integral was evaluated by closing the contour in either the 
upper (z > 0) or lower [z < 0) half-plane. 

Substituting the results into fl2Tj|) we obtain for the average density 



V or f°° r/.A _T_ 



d = — 2L / e- lw -^ |z| *(z)d« = — , (24) 

which coincides with the result of averaging over k with the bond-lengths 
fixed 0,1,0. 
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4 Two-point correlation function 

The two-point correlation function is given by 

R 2 (x) = lim -^R (k,k + = ) , (25) 

where d is the mean density, the limit is taken in such a way that kAL — > oo, 
and we take 

R(k u k 2 ) = (dikjdih))^} (26) 
' f°° -A L r L s e i ^^ {zitanklL ' +Z2t! " nk2L ^ dz \ 

J _ OQ COS 2 fciL r COS 2 fc 2 L s 47T 2 / 



{^} 



with z = (zi, z 2 ). 

In this case, the analogue of ( p0|) is that 



= y {^(z)r'- 1 (z)+,;(t;-l)0 1 (z)0 2 (z)r- 2 (z)}— , (27) 



where 



g( Z ) = — / — — ^ e i(,itan( fcl L)+ 22 tan(fc 2 L)) dL; 

AL J Lo cos 2 k\L cos z k 2 L 

i pLo+AL t 

^( z ) = / e i( Zl t a n(fe 1 L)+^ta I1 (fe 2 i)) dL) (30) 

AL cos 2 fciL 

i pLo+AL j 

(j)J z ) = _ / £ e i(» 1 ta n (*iL)+* a t«i(fc a L)) dL< ( 31 ) 

AL J io cos 2 A; 2 L 

Substituting k± = k, k 2 = k + nx/ (vL ), where x is fixed, and taking the 
limits k — > oo, AL — > (while kAL — > oo), we obtain for the first integral 

y( z ) = J_ |' L0+AL e ^itan(feL)+ 22 tan(feL+^)) rfL 



1 W 2 



7T 



tt/2 



e i( 2l tan^ 2 tan(0+^))^ 5 ^ 



where we have again used L/L fa 1 and, as in the transition from (|2"T|) to 
(|22|), we have approximated / by the integral over one period. We now write 

/ ttx\ tan0 + tan(^f) 1 + /? 2 

tan(0+ — )= = -/? + - 33 

J l-tan0tan(^f) H /5-tan0' v ; 

where (3 = (tan(7rx/v)) _1 oc v/(ttx) (we are interested in the v — > oo limit). 
Performing the change of variables a = tan0 — /?, we arrive at 

i/3{ Zl -z 2 ) roo 2 j 

/( z )-~ / ^ ? /° r - (34) 

Note that /(z) is invariant under the exchange Z\ z 2 and f3 — > —(3, which 
can be verified by the change of variables a = (f3 2 + l)/y in 

To evaluate the integral in (|34}) we differentiate it with respect to Z\ and 
z 2 to get 

df df i e W*-«) /■« __ J22 ^ , _ , /3 2 + ^ da 



e wia -" 2 - 2/3 + a + 



dz\ dz 2 7r J ^ \ a J (a + (3) 2 + 1 



7T 



/CO 2 /-/ 

eizia _, 22 L±iaa = _ e ^i-, 2 ) $ ^ i; ^ ( 35 ) 



where 



ss -- f>™4^ (36) 



./_oo a 



2sign(z 1 )/f(-^ 1 z 2 )J (2v/-(/3 2 + l)^ 2 



Jo(x) being the Bessel function of the first kind and H(x) the Heaviside 
function (characteristic function of the half axis [0, oo)). 
Applying the method of characteristics to the PDE 

|^-|^ = -e^ 1 ^ ) $U,z 2 ), (37) 

OZ\ OZ2 

we obtain the solution 

/(z) = e -l*i+«l _ / 1 e ^(2»-*i-*a)$ Zl + ^ _ y ) dy. (38) 
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Treating the integral for g(z) (see (^)) in a fashion similar to the one 
used to obtain (1341) leads us to 



j2 I"kI2 tan((/!>)+22 tan(0+7rx/«)) 

gM „ ^J^____^^ (39) 



7r ./_oo \ V a 



Comparing this integral to the one in (PB[), and noting that 

l + (l±£. + p)' = ^tlL +p+ t±l +0 ), (40) 

\ a ) a \ a ) 

we have that 

g(z) = LKP 2 + 1) (J- - A) [ e ^--)$(, 1; z 2 )] . (41) 
One can derive a similar expression for the functions 0i(z), 

i/3(2i-2 2 ) /-oo 2 a 

</>i(z) « L / e l2lC ™ » da = L e i/3(21 " 22) — $( Zl , z 2 ), (42) 

7T 7-00 ozi 

and 02 (z), 

2 (z) « L - / + 2 ijdQi 

= -L e^-^)A$( 2l;Z2 ). (43) 

Now we have all the ingredients necessary for evaluating the integral in 
([27|) . Substituting the expression for g(z), (|4T|) , into the first half of the 
integral and integrating it by parts we obtain 

= v{v-l)L% J ^(p* + l)f>-*e H P<.*->*)&. (44) 
11 



Thus 

R 2 (x) 



t 



dz 

4^ 



/ 



v-2 2i/3(z 1 -z 2 ) 



(P 2 + 1)$ 2 



<9$ (9$ 

dz\ dz 2 



(45) 



Now we need to take the limit v — > oo. To do so we write f v 2 (z) = 2 ) ln / 
and rescale /(z) 

1/ ,1 ^7i."-> I / 

/(u//3) = 



ki+"2 
P — 



i r 1 e i (2tf -« 1 -« 3 )^ (j/j Ui + M2 _ (46) 

P Jo 



and hence, to the leading order in 1/(3 = ttx/v, we have 



(v — 2)ln/(u) ~ —7ix | 

= —tixQ, 
where \& is the rescaled function $, 



u 2 - y)dy 

(47; 



2sign(ui)fl"(— «i« 2 ) Jo {2y/-u 1 u 2 ) 



(48) 



and we have taken the limit u — > 00 (/5 — > 00). 

Renormalizing the rest of ([45"D and taking the limit v — > 00 we obtain 



9«i 9«2 



(49) 



The only change when the above calculation is generalized to other boundary 
conditions at the central vertex (i.e. to nonzero values of A -1 in fllq) ) is the 
appearance of a factor e~ x ( Zl+Z2 ) next to every occurrence of dz in the 
above integrals. For A fixed, this factor disappears after rescaling z = u//3 
and taking the limit (3 — > 00. Hence equation ( ]4"9"D is then independent of A. 
In the case when A" 1 = A -1 ?;, the dependence of the spectral statistics on the 
boundary conditions at the central vertex persists. The above expressions 
then coincide with those for those for Seba billiards with a renormalized 



coupling consant, given in 15 



For the derivatives of the function \1/ one has 



dm 

duo 



2 ( J (0)S(ui) + sign(« 1 )iZ"(— uiu 2 



U 2 Jq (2 y /-u 1 u 2 ) 



-J (0)5(m 2 ) + sign(u 1 )iZ"(— u\u 2 
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\J—UyU 2 

uiJ' Q {2^/-uiu 2 



\J—U\U 2 



(50) 
(51) 



therefore, using Jo(0) = 1 and J' Q {x) = —Ji(x), 

d^d^ 2 = ~ 4 i 6 ^ 6 ^ + H {-^ u 2)Jl (2^=^)) • (52) 

Thus 

R 2 {x) = 1 + Je-xQ+Hm-ui) jj2 ( 2> /=^i) + J 2 (2 v ^^)]iJ(-n 1 n 2 )du. 

(53) 

Now we perform the change of variables u 2 i— > — u 2 arriving at the following 
integral representation of the two-point correlation function, 

R 2 (x) = 1 + / e -™^(»)+M(«i-H«) [j2 (2y/w£) + J 2 (2v^)] du. (54) 
Jd 

Here the domain of integration D includes first and third quadrants of the 
MiM2-plane and M(u) is given by 



M(u) = M(ui, u 2 ) = \u x —u 2 \+ / e^-^+^^^Mi-Uz-^dy 

(mi) r (m 2 ) s (r + s-2)! 



«1 + M 2 



izsign^j^ r!s !(r-l)!( S -l)! " (5&j 

r,s=l v ' v ' 

Equation (|54"|) constitutes an exact formula for R 2 (x) for star graphs in 
the limit u — > oo. It is our main result. The point we seek to draw attention 
to is that it is exactly the same as the one obtained in [15] for Seba billiards 
when |?/v(x )| 2 = constant in (0) and A — > oo. We will expand on this 
observation later. First, we consider some of the properties of the two-point 
correlation function and the form factor in more detail. 



5 Expansion for large x 

To derive an expansion of the two point correlation function R 2 (x) for large 
x we notice that since M(— u) = M(u), the integral over the third quadrant 
in (^4j) is equal to the complex conjugate of the integral over second quarter- 
plane, i.e. 

R 2 ( x ) = 1 + 23? / / e-™ M(u)+2i( " 1+U2) J(u)du, (56) 



o 
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where 



J(u) = J 2 (2 v ^i) + J 2 (2^T 2 ) = 



n=0 



-l) n u>£(2n)! 
(n + l)!(ra!)3 



(57) 



Now we can use the expansion of M(u), (0), to expand R2(x) in the powers 
of We substitute itj = ^/(xtt) and obtain 

2i (71 + 72 - 7172) 



R9.X) 



1 + 23? 



£% 2 



■71—72 



1 + 



X7T 

(5 7l72 + 2 7l 2 + 2 7 2 - 5 7l72 2 - 5 7l 2 72 + 2 7l 2 7 2 ) 



x 2 ir 2 



O 



.r 



23? 



2i 



X 2 7T 2 



X 3 7T 3 X 4 7T 4 



(58) 



To compare this to the expansion (|l"4"|) of K{t) we note that if K(t) 
1 + Xlfcli a k rk for r > then, inverting the Fourier transform in @, 



/•oo 

- 1 = 23?lim / (AT(r) - l)e~^ i{x - ie)T dr 



23?V - 
^ \ 2tt 



A'=l 



k+1 



a k k\ 



x 



k+l " 



Applying this to 



K(t) 



1 - 4t + 8r 2 



r r 3 + 0(r 4 ), 



(59) 
(60) 

(61) 



we see that the first few coefficients of the two expansions agree. The proof 
that it is so for all coefficients is given by the following proposition. 



Proposition 1. The asymptotic expansion ( \5dj ) of the two-point correlation 
function and the expansion ( fIZ[ ) of the form factor coincide under the Fourier 
transformation 



PC 

e -nxM(u)+2i(u 1+ u 2 )j^ du = / 

Jo 



(K(t') - 1) e 



- 2niXT 'dr' 



(62) 



Proof. The Fourier transform in (|6~2"|) establishes the correspondence between 
the terms in the asymptotic expansion of 



I* poo 

R 2 {x) = J J e- nxM{u)+2i{ui+U2) J(u)du 



(63) 
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and the terms of the small r expansion of K{r). This correspondence is 

1 r^ 1 



{2mx) k 0-1)!' 



(64) 



Our plan is to modify the integrand in the definition of R 2 (x), getting rid 
of the factor e 2t ( u i+ u 2) J(u), expand the integral in inverse powers of x and 
apply the correspondence rule (|64"D to recover fllPp. 

First of all, as one can verify by direct substitution of the series for 
M(u 1 ,u 2 ), 



® ® ^ ( A/r f ai a2 \\ _ -r+s ( r + A (a 1 /x) r (a 2 /x) 



r,s=0 



(65) 



and 



A ( x m (— ^\) = Y" 2i r+s+1 ( r + s ~ l)K<xi/x) r ((X2/xy 
dx V \ x 1 x J J r\s\(r — l)!(s — 1)! 



gjvgjga ^ f e i(a 1+ a 2 )/,_ ( g 6) 

X \ X / 



Applying (§§, 



(9x 2 

„— TTXM 



(67) 



where <fi = i(a,\ + 0,2) l x an d for simplicity we have omitted the argument 
(ai/x,a 2 /x) of the functions M, J and J\. 
Similarly, using (BH), we have 



d 8 



dcti da 2 

—■kxM 



2 



p \ X ' X j 



(2ixi \ 
4ir 2 J%e 24, (2J Q e <l, aia 2 + iJ^ 'Ja^iai + a 2 )) . 
OL\Ol 2 X v / 
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Noticing the similarity between fl6"T| ) and fl6"8|), we subtract the first from the 
second, with the appropriate factors, to obtain 



4vr 2 



d d 
+ 



x 2 V da\ dct2 



2 



a±a2 dx 2 



= ^ [ ^ + J?] e 2 ^ M , (69) 

where, as before, the argument (cxi/x, a^/x) of M, Jo and J\ has been omit- 
ted. The right hand side of ( |69"D is exactly the integrand of (|56[) if we perform 
the change of variables Uj = ctj/x and, therefore, 



da\da2 
4vr 2 



9 



x 2 V 9«i (9a 2 



5 2 



ai«2 c^ 2 



The first term in the integral can be evaluated as follows, 



(70) 



dotidot2 ( d 







o Air 2 x 2 \dcti da 



V x 7 x / 



da 



4nx 2 



2 r/-^ioo 

ai=0 



da\ rr^ioo 

~ [ w Jo 2 =o ; ' 



2txx 2 



where 

e = 

Since 











dai da 2 



g V x ' x , 



joo _ _2 e ia2/^p-7ra2 
Jai=0 



]oo _ _2„iai/x — ?rai 
«2=0 



we obtain 



da\da2 ( d 







o 47r 2 x 2 \da\ da2 



2vTX 2 TT — i/x 



(71) 



(73) 



(74) 



Now we can expand the result in inverse powers of x and apply the corre- 
spondence rule fl6"4|) . We obtain 

00 '_2 r ) fc + 1 



7TX 7TX — I 



k=0 



fc=0 



(*-l)! 



2 (e" 2r - 1) . (75) 
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Next we need to expand the second part of the integrand in ([701), 

' oo 

2 ™E 



r,s=0 



-e = —e - 7r ( Q i+ Q2 ) exp 



dx 2 dx 2 * \ x r+s+1 r\s\(r + l)!(s + 1)! 



-7r(ai+Q2) . 

dx 2 



?! \ ^ x r+s+1 r!s!(r + l)!(s + 1)! 

j=0 J \r,s=0 \ < J \ < J 



• (76) 



Using the same notation as in 



^ ^r+H^+V + a)! ] _ y (»ai) r+1 (»Q2)' +1 , 



x r+s+1 r\s\(r + l)!(s + 1)! / \^ x 

vr,s=0 v ' v 7 / \r,s=0 



ii,5=0 



where, as before, Fj (R, S) is the jth convolution of F\ (R, S) with itself. Thus 



<9x 2 j! 



E 

i=i 

(R+S+j + l)\x«+SW m ' S) - (?8) 



Finally we integrate against dotida.21 '(47r 2 ! ai«2) to arrive at 

<i«i<i«2 d 2 _„ xM (^i 22.) 
47r 2 a 1 a 2 ^ 2 

^ 4tt 2 j! ^ (R+S + j- l)\(-nr) R + s + 2 ix R+s +i+ 2 A ' J 

j'=l J R,S=0 K J ' 

This is exactly the same as the j sum in ( JlQ| ) with the exception of the extra 
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j = 1 term in the summation above. For j = 1 we have 



2r)' 



(R + sy. » • ^ o ( fl + i) K5+ i)i 



-2r\2 



Js (« +i > ! ilS < s+i > ! r (1 ~ e " 

= 1 - 2e" 2r + e" 4r , (80) 



which, together with the terms 1 and 2(e 2r — 1), gives the correct contribu- 
tion e _4r . □ 



6 Singularities of the form factor 

One can also obtain some information about the singularities of K(t) by 
Fourier transforming the integral representation (0). There is, however, a 
subtle problem associated with this approach. The form factor is by defini- 
tion an even function defined on the real line. What we want to get from 
transforming ( j56|) is an analytic function which coincides with the form factor 
for real r > 0, so as to be able to study its complex singularities. 
As we saw above, 

R 2 (x)= e~ 7TxM{ll)+2iiui+U2) J(u)du= (K(t') - l)e~ 27TiXT '. (81) 

J Jo Jo 

Integrating ( |81~|) against e 2mXT on the real line we obtain 

R 2 {x)e 2nlXT dx = K{t) - 1, r > 0. (82) 



One can check that this leads to the correct power series expansion of^the 
form factor: give x a small negative imaginary part, x h- > x — ie, in R2(x) 
(this is consistent with (|8"ID), substitute in the asymptotic expansion (56), 
and integrate term-by-term. 

We now use R 2 (—x) = R 2 (x) to write 

j e 2niXT R 2 (x)dx = J (e 2 * iXT %(x) + e- 2 ™ XT R 2 (x)) dx. (83) 
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The only factor i? 2 (ic) which depends on x is e 7rxM ( u ) and 

1 



e 27ria;T e -7ra;M(u) ^ 



?r(M(u) - 2ir) : 



(84) 



thus we have for the form factor 
K(t 



7T 



g2j(ui+n 2 ) g 
+ 



J(u)du. (85) 



M(u) - 2zr M(u) + 2ir 
The representation fl8"5]) presents us with a way to find the singularities of 



K(t). These are given by the condition r = M(u s )/(2i) and r = M(u s )/(2z), 
where the point u s is such that 



dM 



dM 



du\ du 2 
The derivative with respect to u 2 is 
dM 



U.s 



0. 



"1 



(86) 



i7) 



where z = ?/ — ui + u 2 and we have assumed that u\ > u 2 > 0. It is obvious 
from the expansion fl5"5j) , however, that the function M(u) is continuously 
differentiable if uiu 2 > and hence that the expression fl87|) is valid for all 
Ui > and u 2 > 0. The integral in ( |8"TD is not easy to analyse and to simplify 
it we reduce our search to the line u 2 = u±, where 



dM 



\u 2 = u x , 



e 2iy J 1 (2y)dy + 2i I e* y J (2y)dy. 



Performing the second integration by parts, 



e 2iy J (2y)dy 

we obtain, after simplification 

dM 



e 2iy U2y) 



2i 



2 
27 



e^J^dy, 



8Ur 



(u 2 = Ul ) = e 2iUl J (2 Ul ) 



(89) 



(90) 



Since §^-(^2 = u i) = §^( u 2 = u i), we see that the zeros of the derivatives 
of M(u) on the line u 2 = u\ are given by the zeros of the Bessel function J . 
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The nearest zero is at u s ~ 1.202. Thus one of the singularities of K{r) lies at 
r s = Af (1.202, 1.202)/(2z) = 0.462 -0.420z. We note that \r s \ = 0.624, which 
coincides with our previous numerical estimate of the radius of convergence 
of the series expansion of K(t) in powers of r around r = 0. This strongly 
suggests that this singularity is the closest to the origin. To this end, we can 
prove the following. 

Proposition 2. Among the singularities arising from stationary points of 
M(ui, u 2 ) along the line u 2 = U\, the singularity at t s = M(1.202, 1.202) / (2i) 
is the nearest to the origin. 

Proof. To show that the statement is true we need to prove that the function 
\M(u,u)\ 2 is a nowhere decreasing function of u. On the line U\ = u 2 = u 
we have 

M(u,u) = / e iy J (y)dy = 2e 2iu u (J (2u) - U^u)) . (91) 
Jo 

Thus \M(x/2, x/2)\ 2 = x 2 (Jq(x) + Jf(x)) and its derivative is, after simpli- 
fication, £\M(x/2,x/2)\ 2 = 2xJ 2 (x) > 0. □ 

It is straightforward to approximate the behaviour of K{r) near these 
singularities. We expand 

M(u) w M(u s ) + --^(u s )(u 1 -u 2 ) 2 + --^(u s )(u 2 -u s ) 2 
d 2 M 

+ t: — ^— (u s )(ui - u s )(u 2 - u s ) 
OU1OU2 

= M{u s ) + a s ({ Ul -u s ) 2 + {u 2 -u s ) 2 ) . (92) 

For the singularity associated with the first Bessel zero, a s ~ 0.385 — 0.349z. 
Then, when r is real, 

K f T) nJ_(r J(u)e 2 ^+^u 

V ; 7ra s J J o ( Ul - u s ) 2 + (u 2 - u s ) 2 + (M(u s ) - 2ir)/a s V ; 

The main contribution to the integral around these singularities is 
K{t) oc -Cln ( 1 - -^-) -C\n(l + JZ= 
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n 



Figure 3: The coefficients of the power series expansion of K{r) normalized 



by p n (crosses), compared to (]95). As expected, the agreement improves as 



n increases. 



where C = J(u s )e 4iUs fa s . Expanding (p4|) into a series around r = we get 



K{t) oc 25? ( cjr P n ^r n ) = 2A^cos(<M + ^)^-r n , (95) 



n=l 



71=1 



where, for the singularity analysed above, A = \J(u s )e AlUs /a s \ ~ 0.519, 
V> = arg ( J(u s )e 4ms /a s ) « -0.737, p = |2z/M(u s )| « 1.602 and = 
arg (2z/M(u s )) ~ 0.737. By Darboux's Principle, the coefficients of the ex- 
pansion (|95|) should comprise the leading contribution to large-order asymp- 
totics of the exact coefficients given by fllHD and (O). To compare them 
we plot the exact coefficients na n / p n against the approximate coefficients 
2Acos((pn + tp). The result is shown in Fig. 0. 
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7 Small x limit of Ro(x 



Returning to (f|9j), one can check that the function defined by (|48D, satisfies 
the equation 



" d 2 


if 9 


- a )] 


2du\du 2 




du 2 J _ 



f^,i{ui-u 2 )^,2\ 

_ 2i(u 1 -u 2 ) 



<9* d^i 



dz\ dz 2 

Substituting it into (|49"D and integrating by parts we obtain 



- * 2 . (96) 



-- I due' 
4 



TTXQ 


r a 2 


if 9 


- 9 )} 




2du\du 2 




du 2 J 



du 



e 2i(ui-n 2 )v]>2 





- 9 ) 


d 2 1 


\dU! 


du 2 J 


2du\du 2 



Now, using the identities 



dui du 2 



d 2 Q 



(e- wxQ ). (97) 



(98) 



.j e i(«i-«a)tt 



29mi(9m 2 

which one can derive using the series expansion of Q(u\, u 2 ) = M(ui, —u 2 ) 



d 2 



we write 








d 




\dU! 


du 2 




= e 


-TTXQ 1 



-11XX 



dQ dQ \ nx d 2 Q (tvx) 2 dQ dQ 



du\ du 2 J 2 du\du 2 2 du\ du 2 



V 2 



2 dui du 2 



Thus we obtain, finally, 

Riix) = - 



^ U c 2i(u 1 -u 2 )-nxQ^ f 2 



7r 2 x 2^_^_ + 3i7rxSS , e i(u 1 -u 2 ) 

OU\ du 2 



(99) 



. (100) 



From ( |100| ) one can see that the two-point correlation function R 2 (x) is linear 

(101) 



in x for small x. The slope was computed in [|i5 

nV3 



R 2 ( x ) = -^- x + 0(x 2 ). 
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8 Discussion 



The derivation presented above provides a proof that two-point spectral cor- 
relations for certain Seba billiards and quantum star graphs are the same, 
in the appropriate limits. This initially surprising fact has its explanation in 
the following observations. First, the dynamics in both systems is centered 
around a single point scatterer; in star graphs it is the central vertex, and 
in Seba billiards the singularity. Furthermore, in between scatterings the 
dynamics is integrable in both cases. 

Second, applying the Mittag-Leffler theorem to the meromorphic function 
tan z, we have that 

00 / 1 1 \ 

tan^= V ; . (102) 

^ \rvK + 7t/2- z nn + n2 v ' 

n=—oo x ' ' 

We can therefore rewrite (|16D in a form similar to ([!]) when |^ n (xo)| 2 = 
constant. It thus becomes less surprising that the two point correlation 
functions of the two systems are the same, because in the limit v — > 00 the 



poles in ( |i5|) have properties similar to those of a Poisson sequence. 

Third, from the mathematical point of view star graphs and Seba billiards 
are similar in that in both cases the scattering centre corresponds quantum 
mechanically to a perturbation of rank one. 

Finally, we remark that our results demonstrate that, at least as regards 
the special case considered here, graphs are able to reproduce features of 
other, experimentally realizable, quantum systems, and also that they pro- 
vide further confirmation that spectral statistics can be computed exactly 
using the trace formula when the periodic orbit statistics are known [|12j. 
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